CP violation in Bd,s l^l' in the model III 2HDM 
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We have calculated the Wilson coefficients Cio,Cq- (i=l,2) in the AIS renormalization scheme 
in the model III 2HDM. Using the obtained Wilson coefficients, we have analyzed the CP violation 
in decays Bq (q=d,s) in the model. The CP asymmetry, Acp, depends on the parameters 

of models and Acp in Bd l^l^ can be as large as 40% and 35% for I = r and I — fi respectively. 
It can reach 4% for decays. Because in SM CP violation is smaller than or equal to O(10~'^) 
which is unobservably small, an observation of CP asymmetry in the decays Bg — > l'^l~(q — d,s) 
would unambiguously signal the existence of new physics. 

PACS numbers: 11.30.E, 13.20.H, 12.60.F, 12.60.J 



I. INTRODUCTION 



OO ' The flavor changing neutral current process, Bd,s — > (l=/i,r), has attracted a lot of attention since it is very 
sensitive to the structure of SM and potential new physics beyond SM and was shown to be powerful to shed light on 
' the existence of new physics before possible new particles are produced at colHders [E S S S H • For example, in a 
: very large region ofparameter space supersymmetric (SUSY) contributions were shown to be easy to overwhelm the 
' SM contributionP,|2,|3,0)l3 and even reach, e.g., for l=/x, the experimental upper boundQ 

B^(Bd ^ /i+Ai") < 6.8 X 10"^ (CL = 90%) 

Br{Bs fi+fj.-) < 2.0 X 10"^ {CL = 90%). (1) 

In other words measuring the branching ratio of Bd,s — l^l" can give stringent constraints on the parameter space of 
I new models beyond SM, especially for that of the minimal supersymmetric standard model (MSSM) because of the 
' tan'^/3 dependence of SUSY contributions in some large tan/3 regions of the parameter spaceQ,|3,l3jl3- 1^ the model 
. II two Higgs doublet model (2IIDM), the branching ratio of Bs — > A*~'_M~ can also reach 2 x 10~^, which is within 
ii' good reach at Tevatron Run II, if tan/3 is large enough (say, > 60 ) 0,0. Comparing with hadronic decays of B 
Oh, mesons, this process is very clean and the only nonperturbative quantity involved is the decay constant that can be 
^ ' calculated by using the lattice gauge theory, QCD sum rules etc. 

An observable without hadronic uncertainty at all is the CP asymmetry in the process Bd^s — * l^l^ since the 
J> ■ common uncertain decay constant, which is the only source of hadronic uncertainty for the process, cancels out. CP 
' violation in the 6-system has been established from measurements of time-dependent asymmetries in _B — > J/'^K 
. decays 0,^3- It is urgent and important to study CP violation in more processes, including the process Bd.s l^l^ 
d ' which might be measured in the near future. Obviously for the process Bd.s l^l^ there are no direct CP violations 
since there are no strong phases in the decay amplitude. But it is well known that CP violating effects can survive 
after taking into account the mixing of the neutral mesons, and in the absence of the strong phases. Recently, 
it has been shown that CP violation in Bd,s — > l^l^ (1=Mj ''') is also an interesting observable for searching for new 
physics jO, the present paper we study it in the model III 2HDM 0, It is well-known that in the 

model III 2HDM the couplings involving Higgs bosons and fermions can have complex phases, which can induce CP 
violation effects, even in the simplest case in which all tree-level FCNC couplings are negligible. So it is expected 
that a significant CP asymmetry in Bd,s 0-=IJ-,t) should exist due to the effects of such extra phases in the 

model. Another motivation to consider the model III is that one starts in the minimal supersymmetric standard 
model (MSSM) with a Higgs sector of type II 2HDM and, after integration of squarks and gluinos, one ends up with 
an unconstrained model HI 2HDM 15]. The Wilson coefficients of operators relevant to b ^ sl'^l~ in the effective 
Hamiltonian have been calculated in the model III 2HDM using the on-shell renormalization prescription . Some 
Feynman diagrams containing the coupling of a neutral Higgs boson to the charged Higgs and Goldstone bosons were 
missed in the calculations in ,16.] . We calculate them, including all the contributions, based on the AI S renormalization 
scheme. We then analyze the CP asymmetry in Bd^s l^l^ ''') using the Wilson coefficients obtained by us. It 

is shown that the CP asymmetry can be as large as 40% for i?° and 4% for i?° in the reasonable region of parameters 
in the model. 

The organization of the paper is as follows. In Section II we describe the model HI 2HDM briefly. In section III 
we give the effective Hamiltonian responsible for b sl~^l~ in the model. We present the formula for CP asymmetry 
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in Bd.s — > in Section IV. The Section V is devoted to numerical results. In Section VI we draw conclusions and 
discussions. Finally we give contributions to Wilson coefficients Cio,Cq. from individual diagrams in the model III 
2HDM in the appendix. 



II. THE MODEL III TWO HIGGS DOUBLET MODEL 



In general one can have a Yukawa Lagrangian of the form 



(2) 



where Hi, for i=l,2, are the two scalar doublets of a 2HDM, while rf'^ and ^^'^ are the non-diagonal matrices of 
the Yukawa couplings. 



For convenience we can choose to express Hi and H2 in a suitable basis such that only the Tj^-^ couplings generate 



the fermion masses, i.e. such that 



{Hi) = ( ^ ) . (^2) = 



(3) 



The two doublets in the basis are of the form 



Hi 



1 

71 



x/2G+ 



i^2 



_ V2H+ 
V2 V <^2 + 



(4) 



where G°'^ are the Goldstone bosons, H^ and are the physical charged-Higgs boson and CP-odd neutral Higgs 
boson respectively. The advantage of using the basis is that the first doublet Hi corresponds to the scalar doublet of 
the SM while the new Higgs fields arise from the second doublet H2- So do not have couplings to the gauge bosons 
of the form (jP^ZZ or (^IW+W' . 

In eq. Q 0°, 02 are not the neutral mass eigenstates but linear combinations of the CP-even neutral Higgs boson 
mass eigenstates, H^ and 



i/" = (i? cos a + 6^ sin a 



(5) 



i)\ sin a -I- (/)2 cos a 



where a is the mixing angle, such that for a — Q, (0]', coincide with the mass eigenstates. 

After diagonalizing the mass matrix of the quark fields, the flavor changing (FC) part of the Yukawa Lagrangian 
becomes 



Cy.FC = iYJQ^,LH2U,.R + £.?AxH2D,,R + h.c. (6) 

where (5i,L, tj,_R, and Dj^R now denote the quark fields which are in mass eigenstates and 

In eq. V^'^ are the rotation matrices acting on the up- and down-type quarks, with left or right chirality 

respectively, so that Vckm = {VlV^l' usual Cabibbo-Kobayashi-Maskawa (CKM) matrix. Feynman rules of 
Yukawa couplings follows from eq. ^ and can be found in, e.g., ref. Il8j|. The flavor changing neutral current 
(FCNC) couplings are given by the matrices ^^'^ and the charged FC couplings are given by 

enlarged = ' VcKAI 

Charged = VcKM ' ■ (8) 

Because the definition of the ^fj^ couplings is arbitrary, we can take the rotated couplings as the original ones and 

shall write ^^'^ in stead of ^'^'^ hereafter. It is worth to note that the peculiar form of charged FC couplings, eq. 

is of an important distinction from popular model I and II and has significant phenomenological effects different 
from those in models I and II . 
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In order to proceed we assume the Cheng-Sher ansatz 



- ~ 

which ensures that the FCNC within the first two generations are naturally suppressed by small quark masses. In the 
ansatz the residual degree of arbitrariness of the FC couplings is expressed through the Xij parameters which are of 
order one and need to be constrained by the available experiments. In the paper we choose ^^'^ to be diagonal for 
the sake of simplicity so that besides Higgs boson masses only Xu and Xbb in the quark sector and A;; in the lepton 
sector are the new parameters which enter into the Wilson coefficients relevant to the rare leptonic B decays. 

The Yukawa Lagrangian Q has no discrete symmetry. In general, the Higgs potential in 2IIDM which has no 
discrete symmetry may be CP conserved or CP violated [iJi 123 • In order to concentrate on the effects of FC Yukawa 
couplings we assume the potential is CP conserved. Therefore, the new source of CP violation is only from the Yukawa 
couplings. In a general basis (i=l,2) the general CP conserved potential is given by: 



F($i,$2) = Ai($I$i-W?)2+A2($^$2-«2)' 

+X3[{^\'S>i^vl) + {^l<i>2-vl)]^ 

+ A4[($I$l)($t$2) - ($t$2)(*^$l)] 

+A5[Re($l$2) - «i«2]' + A6[Im($l$2)]' 
+A7($I$i - v?)(i?e($t$2) - V1V2) 

+A8($J$2 - v|)(i?e($l$2) - V1V2) . (10) 
in which all the coupling constants Xi are real due to hermiticity. The minimum of the potential is at 

< $1 f \ < $2 >= f M , (11) 



thus breaks SU{2) x U{1) down to U{1)em- For the sake of simplicity and to decrease the number of the parameters 
we assume Ai = A2 and A7 = Ag = 2 A3 + A5/2 in the potential H10() . Then we have seven parameters altogether in the 
Higgs sector. One of them, v — y^vf + w|, is fixed by the W boson mass, = ^g^v^ . The others can be expressed 
in terms of mixing angles a' and (3 and Higgs masses mfj± , m^o , m^/o , 77*^,0 . From the potential, it is straightforward 
to derive 

9A0H-G+ — 7; > U^j 



9H0H-G+ 



9hOH-G+ 



2mw 
igsui{P — Q;')(m^o 



''H±) 



2mw 

"h° ~ "''H±) 



-igcoa{(3 — a')(m?o — m 



2m 



w 



which have the same forms as those in a general model II 2HDM |J, |21| . 

The specific basis H2), eq. Q), can be obtained by a unitary transformation 



H = S<^, S=[ ) > (13) 



— sin/3 cos/3 

where H = {Hi,H2)'^, <f> = (<I>i, $2)^, and tan/3 = V2/V1. 

III. EFFECTIVE HAMILTONIAN FOR b s{d)l+l- 

The effective Hamiltonian describing the flavor changing processes b s{d)l^l^ can be defined as 

10 10 

H^ff = -—^XtiY, aii^mf^) + cq, (m) 

^ i=l i=l 



where Xt = VtbV^*^, 0[s{i = 1, • • • , 10) are the same as those given in the ref.|22,|2^, and Q[s come from exchanging 
neutral Higgs bosons and have been given in refs. pll2^. 
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The explicit expressions of relevant operators in the Model III 2HDM for the process Bq — > (q=d,s) are 

Oio = ^(qti.btKh.Toi) 

Qi = ^imm 

We calculate their Wilson coefficients at one loop level in the MS renormalization scheme and results are given in 
Appendix. Before proceeding, a remark is in place. In calculations of the Wilson coefficients Cq. the neutral Higgs 
boson h° ^°)-penguin diagrams with the charged Higgs and Goldstone bosons in the loop, which is missed in 

ref. should be included. In the Appendix we write the contributions of these diagrams to Cq^ in terms of the 
general coupling constants gA°H-G+j 9h°h-g+ gh°H-G+- However, in numerical calculations we shall confine 
ourself to the specific case of the couplings H12(l . 

The QCD corrections to coefficients Ci and Cq. can be incorporated in the standard way by using the renormaliza- 
tion group equations. Qi{i = 1, • • • , 10) does not mix with Oiq so that the evolution of Cio remains unchanged and 
is given by |22| 

Cio(TOfc) = Cio{mw)- (15) 

It is obvious that operators Oi{i — 1, • • • , 10) and Qi{i = 3, • • • , 10) do not mix into Qi and Q2 and also there is 
no mixing between Qi and Q2- Therefore, the evolution of Cq-^,C'q^ is controlled by the anomalous dimensions of 
Qi,Q2 respectively |2J|, i.e., 

= 1.24 Cg^imw), « = 1,2, (16) 
where 7q = —4 is the one loop anomalous dimension of qLbR, 77 = ~ 1-72, and /So = 11 — (2/3)7i/ = 23/3. 



IV. CP ASYMMETRY IN Ba.s ri~ 

We need to know what kind of CP violating observables can be defined in the process Bq — > for q=d, s. At 
first, direct CP violation is absent in this process due to absence of strong phases. T-odd projection of polarization 
is a kind of useful tool to probe the CP violating effects, for example, \n B ^ Xs l'^ l~ [26, 27, 28j. However for the 
process we are discussing here, we have actually only one independent momentum and one independent spin which 
can be chosen as those of l~, so no T-odd projections can be defined. Unlike the case generally discussed for hadronic 
final states, for example, that in Ref. the detected final states of Z+ and l~ of this process in experiments are 
basically two asymptotic energy-momentum eigenstates which are not CP eigenstates. Considering for instance B^ 
decays to in the rest frame of -B", due to the energy-momentum conservation we denote the four-momenta of l~ 
and Z+ as p = {E,p) and p — {E,—p). Then the angular momentum conservation tells us that and l^l]^ final 
states are forbidden. Hence we are left with a pair of CP conjugated final states, and and the couple of 

the corresponding CP conjugated process. Therefore, we may define the time dependent CP asymmetries as 



3upi( 

m 



^CP\t) - TTTT) 77s , T./nn 77 7+7^ y^') 



nB',,yAt)^ipR)+nBXsit)^itii) 

Two corresponding time integrated CP asymmetries are 

"^^^ jrdtnBXAt)^h) + ^rdtnBi,yAt)^h) ' ^ ^ 

Where fi 2 = r^l r with Il{r) being the helicity eigenstate of eigenvalue — 1(+1), / is the CP conjugated state of 
/• 
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The time evolutions of the initial pure and states are given bv|30| 

\B°,,yM)=9+it)\B') + ^9-itW), 
P 



with g± (t) given by 



\B°Hysit))^-^g-{t)\B')+g+it)\B°). (20) 



g+{t) = exp{--rt - imt)cos(—^t), 

1 Am 

g-(t) = exp(~—Tt — imtjisin[—^t) (21) 



The absence of strong phases implies 



\Af\ = |%|, \Af\ = \Af\ (22) 
where Af{Af) =< f\Heff\B°{B°) >■ And the CPT invariance leads to 



Af ( Aj 



Af \Aj 



(23) 



In the above definition of CP asymmetry we need to separate the final state from V^l]^ in order to measure CP 
asymmetry. For l=r, the polarization analysis is straightforward. However, detecting tau's is difficult experimentally. 
For l=/i, in principle one can separate the final state /^J/i^ from by measuring the energy spectra of the 

electron from muon decay A /i^ will decay to an energetic cl, which must go forward to carry the muon spin, 
and a less energetic pair of neutrino and antineutrino because the electron is always left-handed^nofootinbib and the 
energy-momentum and angular momentum are conserved. Due to the same reason, for /i/j, the relative energies of 
electron and a pair of neutrino and antineutrino are roughly reversed. Therefore, the energy spectra of the electron 
from the muon decay is a powerful spin analyzer. However, in practice muons never decay in a 47r detector because 
the lifetime of a muon is long ( ct=659 m). As pointed out in ref.|llj. a possible way to measure a polarized muon 
decay is to build special detectors which can make muons lose its energy but keep polarization so that the polarized 
muon decays can be measured. 

In order to make measurements accessible, as proposed in ref. [l^ . one can define the CP violating observable as 

D 

AcP = ^, 



dt J2 nB^nysit) ^ f^) 
■^0 ^=l,2 

~ dt r(5X.(0 - fi). 

5 = / dt Y nB^sit) ^ /o 

■^0 ,=1,2 



dt J2 r(5X.(0 ^ /.) (24) 

1=1,2 

Different from A]~,p (i=l,2), it is accessible experimental to measure such defined CP asymmetry Acp. So in the 
following we shall concentrate on it. 

Using the effective Hamiltonian, we obtain by a straightforward calculation ^ 



% ^ At Cqiy/l - 4mf + (Cq2 + 2miCio) 



^ In the present case it is quite a good approximation to ignore the mass of electron. 

^ We have neglected the contributions, which is smaller than or equal to 10"'^ of the leading term, from the penguin diagrams with c and 
u quarks in the loop. It is true for both B^i and Bs decays |3^ . Therefore, although there are weak phases from the c or u quark in the 
loop, in particular, for B^, the effect on the decay phase induced by them is negligibly small. 
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where At = Vt^VJ^ or VtbV^*, mi = mi/mB« and C^'s are the Wilson coefficients at the rrn, scale. Because Cg/s are 
proportional to m; and Cio is independent of m; it follows from eq. (|25|l that the CP asymmetry in g —> I'^l^ is 
independent of the mass of the lepton in the approximation Of neglecting Arhf. That is, it is the same for 1 = electron, 
muon. 

In SM, one hasjH^ 

q__Mr^__K .... 
Wh,\- At' ^'"^ 

up to the correction smaller than or equal to order of lO^'^, Cio is real, Cq^ = 0, and Cq^ is negligibly small. So it 
follows from Eas. H25() . (|26|l that there is no CP violation in SM. ^ If one includes the correction smaller than order of 
10^^ to = one will have CP violation of order of lO^'^ for i?° and 10~^ for B'^ which are unobservably small. 
In the approximation | ^ | = 1 the time integrated CP asymmetry is 



where Xq = {q = d,s for B^ and B^ respectively) , and ^ is 

Cqi v/1 - 4.mf + {Cq2 + 2miC 



^CP = - , i^iow-, , V?. ' 1 = d,s, (27) 



CQiyT^W-(CQ2+2™/Qo) 



(28) 



As expected, it is nonzero in the presence of CP violating phases. In deriving Eq. (|27|l we have used eq. GHJ which 
is the result for B*^-B^ mixing in the SM. B^-B'^ mixing in the model II 2IIDM has been examined in ref. In the 
model III there are new CP violating phases which might affect the mixing. However, for the values of parameters 
which we assumed ( see below, ea. (|5n)l ') the correction to the SM value of q/p is below 20% and consequently we can 
still use eq. (I27|l as an approximation. Because Xg is larger than 19.0(90% CL) and Xd is just about 0.76 [a| Acp in 
Bs decays is much smaller than that in Bd decays, as can be seen from eq. (|?7|l. 



V. NUMERICAL RESULTS 



In numerical calculations, the following values of parameters are assumed: 

ruAo = 120Gey, ruho = llbGeV, m^o = IGOGeV, mH± = 200Gey . (29) 

To study the dependence of Acp on masses of Higgs bosons, we also include below some results obtained with the 
doubled values of eq. H29|l for Higgs boson masses. 

The plots of |CqJ and \CqJ for r versus |A6f,| for IAttI — 5 and fixed |Att| are shown in Fig. 1 and Fig. 2 
respectively. One can see from the figures that |CqJ increases when |Abf,| increases and the increase is faster for a 
larger |Att| than that for a small |Att|. Comparing with the model II 2HDM, |CqJ can reach the values larger than 
those in the model II if the value of |Af,b|, taken to be >> 1, is equal to that of tan/3 in the model II and the parameter 
|Att| is larger than l/|A;,b|. We calculate the constraint on \Xbb\ and |Att| for jA^^^j = 50 due to the experimental upper 
bound of Br{Bs — > /i^/i~) and the result is shown in Fig. 3 where the horizontal line represents the upper limit of 
|Att| which comes from the experimental constraints oi B — B mixing , r(6 ^ 57), r(6 crDr), Po, Rb and electric 
dipole moments (EDMs) of the electron and neutron p^ . The region under the curve (and the horizontal line) is 
allowed by the experimental bound of Br(i?s fi^ fi~). We also calculate the corresponding constraint by including 
the phases of the two parameters. The result is that the constraint is not sensitive to the phase of Xbt or Att, as 
expected. 



^ Note that the phase convention between and is fixed as CPI-B" >= — |B" > when deriving eqs. I25i . 1261 . 

One can checlt by combining Eqs. I26i and I25i that all freedoms of phase conventions are cancelled out completely in --j^, including 

P fi 

the one between B'^ and B^. 

^ According to the box diagram calculation in SM, the deviation of \q/p\ from 1 is 10~^(10~^) for B^( B.,) ^^ . So 10~^ is a conservative 
estimate. 
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In numerically calculating Acp, the following values of parameters are assumed: 

|Af,b|=50, |Au|= 0.02, (30) 

|A,,| = 5,50, |A^^|=5,50, Orr^Tr/A, e^^^n/A, (31) 

and 

ett+0bb = 7T/2,TT/3 (32) 

has been set for simplicity. The values of Higgs boson masses and Xbb.tt, i.e., eqs. H29|l . (jSOJ, and (|32|l have been set 
to satisfy the experimental constraints oi B — B mixing , T{b — *■ 57), T{b — > ctDt), pq, Rb and electric dipolc moments 
(EDMs) of the electron and neutron ^^l- The value of IAttL eg . H31|l . has been set to satisfy the constraint from 
Z l^l^ which has been analyzed and the upper limit is 56 The Figs. 4 and 5 are devoted to Acp versus Ou 
for B^ T^T^ and i?° M"*"/^" respectively. Acp depends on 6^ significantly and can reach 40% and 35% for / = r 
and / — ji respectively. As for the B^ decays, Acp is much smaller than that in B^ decays, e.g., it can only reach 
4% for T^T^ , as expected. In order to see the effects of different parameters we consider two cases (a) and (b) 

which correspond [At-^I (or |Ap^|) =5 and 50 respectively and plot two curves in each case in Figs. 4,5. In Fig. 4 the 
solid (dotted) one corresponds the set of Higgs boson masses same as eq. (29) (doubled values of those in eq. (29)). 
One can see from Fig. 4a that the maximum of Acp significantly decreases when the masses of Higgs bosons increase. 
However, Fig. 4b shows that the change is small when the masses of Higgs bosons are doubled. The reason is that in 
the case (a) \Cq- \ (i=l,2) for low Higgs masses, eq.(29), is the same order of magnitude as rhrCio and Cio is almost 
real, when Higgs boson masses increase Cg. decreases but Cio keeps unchanged so that the CP asymmetry decreases, 
as shown in Fig. 4a. In the case (b), Cq- is much larger than mrCio. Therefore, when Higgs boson masses increase 
although Cq. is decreased ^ is almost unchanged since the numerator and dominator are almost simultaneously scaled 
(see, eq. (|28|l ). The small change of Aqp is due to the increase of m^±. In order to see the effect of the phase 
Ou + Obb on Acp we plot two curves in Fig. 5 where the solid (dotted) one corresponds 9u + Obb = 7i'/2 {tt/3). One can 
see from the figure that the curve moves toward the direction opposite to that of the transverse axis when 0tt + Obb 
decreases. 

With the branching ratios 



l-imf)Ci,^+{CQ,+2miCiof , (33) 



where tb^ is the Bq lifetime, we calculate the events Nq needed for observing Acp at la in the areas of parameter 
space in which Acp and the branching ratios both have large values and all experimental constraints are satisfied. 

|io_(l 

the CP asymmetry in i? ^ with good accuracy. For l=r, the events Nq are order of 10^ and 10^ for B^ and 

5° respectively. Assuming a total of 5 x 10^(10^) B^Bd [BsBs) decays, one can expect to observe ~ 100 identified 
Bq — !• r+T~ events, permitting a test of the predicted CP asymmetry with good accuracy. 



For l=/i, they are order of 10^ and 10^ for B^ and i?° respectively. Therefore, 10 ^° ( 10^^) Bd [Bg) per year, which 
is in the designed range in the future B factors with 10^- 10^^ B hadrons per year [23, is needed in order to observe 



VI. CONCLUSIONS 



In summary, we have calculated the Wilson coefficients Cio, Cq. (i=l,2) in the MS renormalization scheme in the 
model HI 2HDM. Comparing with the model II 2HDM, |Cq. | can reach the values larger than those in the model II 
when the value of |Afcb|, taken to be >> 1, is equal to that of tan/3 in the model II and the parameter |A(t| is larger 
than l/|Afcb|. It is shown that there is a constraint on Abb and Xu due to the experimental limit of Br(i3s — > ji^ ^^). 

We have analyzed the CP violation in decays B^ (q=d,s). The CP asymmetry depends on the parameters 

of models, in particular, the phase Ou significantly. Aqp in Bd can be as large as 40% and 35% for Z = t 

and / = respectively. It can reach 4% for B^ decays. Because in SM CP violation is smaller than or equal to 
O(10~'^) which is unobservably small, an observation of CP asymmetry in the decays B^ l^l^{q — d,s) would 
unambiguously signal the existence of new physics. 
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Appendix 

By computing the self-energy type, Higgs-penguin and box diagrams, Cq^ and Cq^ for / = t with the superscript 
denoting the type of a diagram or the contribution of exchanging a specific Higgs boson are extracted out and given 
below 



XtmbnirlmXrTi 



1 pl-X 

dx / dy 
Jo 



ytXui'2-x) ^ ytXttjl + x) 
Si{yt,xt) Si{xt,yt) 



+ \Xtt\'^Xbb\l + 2\n 



S2{yt)mjj- _^ ytiXtt - 2iImXttx) ytXu 



+ l + 21n 



S2{xt)7 



'W' \ + ^ ( l + 21n 



Xt 



m%_Si{yt,xt) 



XtlJ." 



_^ '^ytgAOH-G+ f ^bb-X^fX ^ Xtt{l-x)\ 2[Xlf + iImXtt{2 + xt)x] 



+ ^ttXlh 



gmw \Si{yt,xt) Si{xt,yt) 

2 pi 

+ Xbb 



S2{Xt) 



In-^^ / dx\nS4,{yt) + Xbb In -^^g / dx\nSi{xt) > 

"^Jf- Jo J L Jo J J 



(34) 



^0 XtmbTrirReXT 

^Q2 = 



4m^o sin^ 



1 pI-x 

dx I dy 





ytX;t{2-x) ytXtt{l+x) 



+ |A„|'A66 l + 21n 



w {Jo Jo I Si{yt,xt) Si{xt,yt) 

S2{yt)m]j. yt{Xtt - 2iImXttx) ytXt 



+ 



Kt I l + 21n 



S2{xt)m 



^(l + 21n 

Xt 



S2{yt)Xbb S2{yt)X*t 

mjj^Si {yt,xt) 



+ 



/ Xt \ XtlJ.^ 

'^yt9A"H-G+ f ^bb - KtX ^ Xtt{l-x)\ 2[Xlf + iImXtt{2 + Xt)x] 



+ ^ttXbb 



gmw \Si{yt,Xt) Si{xt,yt) 

2 rl 

2 ' 



S2{xt) 



In 



m 



H- 



dx\nS4:{yt) 







+ Xbb 



In- 



m 



w JO 



dx In S4 {xt ) 



(35) 



Xtmbmr{ReXrTSa + Ca) 



4m? 



dx 
Jo 



dy 



2{Sa{X*t - 2ReXttx) + (7^(1 - 2a:)) 

S2{xt) 



+ 



SaytXlt{x - 2) 2xt{ReXttSa + Ca)x ytXttSa{l + x) Ca{4:X-3xt) 



Si{yt,xt) 
(A*t^„ + a) ( l + 21n 



^('l + 21n'"^ 



S2{xt) 

rnhS2{x^ 



Si{xt,yt) 



+ 



xtSsixt) 



SaXbb /^j^^21n™^"'^^*'^*'^*' 



Xt 



XtfJ.^ 



yt{So,{Xtt - 2ReXttx) + ga(l - 2x)) 
S2{yt) 
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/ dx In 54(2/4) - In— 5- 
Jo 
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Jo i^w J J 
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c: 
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dx 
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dx 
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1 



Siiyt,xt) 
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Adding all the contributions to Cg^ and Cq^ given above respectively, we get the Wilson coefficients Cq- (i=l,2). 
In numerical calculations we use the couplings among the neutral and charged Higgs bosons and charged Goldstone 
bosons given in eq. H12I) . 

The Wilson coefficient Cio is 



Cio — 



Cnew 
10 



10 + '-'lO ; 

2j/t(3-4sin^ 9w) 
S2{yt) 

(3- 2sin2 Ow) 



l-x 



dy 







6(2 sin^ Ow ~ l)ln 



— » sm tiw 

1 



1 + ln 



m]j_S2{yt) 



dx{x — 1) lnS'4(2/f) + ■ 







"H- 



(42) 



In above equations, the definitions of the functions Si are 



SiipA) = p(t){l ~ X - y) + x<j) + yp 

S2{p) = {l-x-y) + {x + y)p 

Saip) = {l~x~y)p+{x + y) 

Siip) = 1-x + xp, (43) 



and 



xt = = , (44) 

Cq = COS a, Sa = sin a, (45) 
with a being the mixing angle of the CP-even neutral Higgs bosons. 



[1] C.-S. Huang and Q.-S. Yan, Phys. Lett. B442 (1998) 209; C.-S. Huang, W. Liao and Q.-S. Yan, Phys. Rev. D59 (1999) 
011701. 

[2] S.R. Choudhury and N. Gaur, Phys. Lett. B451 (1999) 86. 
[3] K.S. Babu and C. Kolda, Phys. Rev. Lett. 84 (2000) 228. 

[4] C.-S. Huang et al., Phys. Rev. D63 (2001) 114021; ibid. 64 (2001) 059902(E); P.H. Chankowski, L. Slawianowska, P.R. 

D63 (2001) 054012; C. Bobeth et al., Phys. Rev. D64 (2001) 074014; C. Bobeth et al., Phys. Rev. D66 (2002) 074021; 

G. Isidori, A. Retico, JEEP 11 (2001) 001. 
[5] A. Dedes, H. K. Dreiner and U. Nierste, Phys. Rev. Lett. 87 (2001) 251804; R. Arnowitt et al., Phys. Lett. B538 (2002) 121- 

129; Z. Xiong, J. M. Yang, Nucl. Phys. B628 (2002) 193; S. Baek, P. Ko and W. Y. Song, hep-ph/0205259 H. Dreiner, U. 

Nierste and P. Richardson, hep-ph/0207026 J. K. Mizukoshi, X. Tata and Y. Wang, hep-ph/0208078 For a recent review, 

see, e.g., Chao-Shang Huang, hep-ph/0210314 a talk presented at 5th International Conference on Hyperons, Charm and 

Beauty Hadrons, Vancouver, Canada, June 25-29, 2002. 



11 



[6] PDG, Review of Particle Data 2002. 

[7] C. Hamzaoui, M. Pospelov and M. Toharia, Pliys. Rev. D59 (1999) 095005. 

[8] H.E. Logan and U. Nierste, Nucl.Phys.B586 (2000) 39-55. 

[9] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 87, 091801 (2001). 
[10] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 87, 091802 (2001). 
[11] C.-S. Huang, Liao W., Phys. Lett. B525 (2002) 107. 
[12] C.-S. Huang, Liao W., Phys. Lett. B538 (2002) 301. 

[13] T.P. Cheng and M. Sher, Phys. Rev. D35 (1987) 3484; D44 (1991) 1461. 

[14] W.S. Hou, Phys. Lett. B296 (1992) 179; A. Antaramian, L.J. HaU, and A. Rasin, Phys. Rev. Lett. 69, 1871 (1992); L.J. 

Hall and S. Weinberg, Phys. Rev. D 48, R979 (1993); M.J. Savage, Phys. Lett. B266, 135 (1991); L. Wolfenstein and 

Y.L. Wu, Phys. Rev. Lett, bf 73 (1994) 2809. 
[15] A. Dobado, M.J. Herrero and D. Temes, Phys.Rev. D65 (2002) 075023; A. Pilaftsis, C.E.M. Wagner, Nucl.Phys. B553 

(1999) 3; M. Carena, J. Ellis, A. Pilaftsis, C.E.M. Wagner, Nucl.Phys. B586 (2000) 92; A. Dedes and A. Pilaftsis, 

hep-ph/0209306 

[16] T.M. Alievluialvi. Savci, Phys.Lett. B481 (2000) 275-286; E.G. Rtan and G. Turan, Phys.Rev. D63 (2001) 115007. 

[17] D. Atwood, L. Reina and A. Soni, Phys.Rev. D55 (1997) 3156. 

[18] D. Bowser-Chao, K. Cheung, and W.-Y. Keung, Phys. Rev. D59 (1999) 115006. 

[19] Z.J. Xiao, C.S. Li, and K.T. Chao, Phys. Lett. B473 (2000) 148, Phys. Rev. D62 (2000) 094008; T.M. Aliev, M. Savci, 
hep-ph/0003188 

[20] F.J. Botella and J.P. Silva, Phys. Rev. D51 (1995) 3870-3875; H. Georgi, Hadronic Jour. 1 (1978) 155. 

[21] J. F. Gunion, H. E. Haber, G. Kane and S. Dawson, The Higgs Hunter's Guide (Addison- Wesley, Reading, MA, 1990); 

errata hep-ph/9302272 
[22] B. Grinstein, M.J. Savage and M.B. Wise, Nucl. Phys. B319 (1989) 271. 
[23] G. Buchalla, A.J. Buras and M.E. Lauthenbache, Rev. Mod. Phys. 68 (1996) 1125. 
[24] Y.-B. Dai, C.-S. Huang, and H.-W. Huang, Phys. Lett. B390 (1997) 257. 
[25] C.-S.Huang, Commun. Theor. Phys. 2 (1983) 1265. 
[26] F. Kriiger and L.M. Sehgal, Phys. Lett. B380 (1996) 199. 

[27] C-S Huang and LIAO Wei, Phys. Rev. D61 (2000) 116002; Phys. Rev. D62 (2000) 016008. 
[28] C.-S. Huang, Z.-H. Zhu, Phys. Rev. D61 (2000) 015011; Erratum-ibid. D61 (2000) 119903. 
[29] M. Gronau, Phys. Rev. Lett. 63(1989) 1451. 

[30] LI. Bigi et al., in "CP violation", ed. C. Jarlskog (World Scientific, Singapore, 1989). 

[31] A similar analysis was presented for seperating tL from tR. See, C.R. Schmidt and M.E. Peskin, Phys. Rev. Lett. 69 (1992) 
410. 

[32] Y. Nir, hep-ph/9911321 (IASSNS-HEP-99-96). 

[33] A.J. Buras, in: Probing the Standard Model of Particle Interactions, Edited by F. David and R. Gupta, Elsevier Science, 

Amsterdam, 1999, hep-ph/9806471 
[34] A.J. Buras, P.H. Chankowski, J. Rosiek, and L. Slawianowska, Nucl. Phys. B619 (2001) 434. 
[35] E.G. Iltan, Phys. Rev. D65 (2002) 036003. 



12 



30 - 



25 - 



20 - 



15 - 




FIG. 1: |CqJ versus \\bi,\. 




FIG. 2: |CqJ versus |A6i,|. 
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FIG. 4: (a) Acp for — > t~^t~ versus the CP violating phase 9tt, for \Xtt\ = 5. The solid line stands for the masses of Higgs 
bosons same as eq. (29), the dotted line stands for the masses of Higgs bosons same as doubled those in eq. (29). (b)The same 
as (a), except for lArrl = 50. 
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FIG. 5: (a) Acp for versus the CP violating phase 6tt, for |Af,^| = 5. The solid line stands for Ou + Obb = 7i'/2 , 

the dotted line for 9tt + Obb = 7r/3. (b)The same as (a) except for |A|uy^| = 50. 



